Nonequilibrium domain formation by pressure fluctuations 
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Fluctuations in thermal many-particle systems reflect fundamental dynami- 
cal processes in both equilibrium and nonequilibrium (NEQ) physics. In NEQ 
systems Jl| fluctuations are important in a variety of contexts ranging from pat- 
tern formation to molecular motors |jj-l7|. Here, we address the question 
if and how fluctuations may be employed to characterize and control pattern 
formation in NEQ nanoscopic systems. We report computer simulations of a 
liquid crystal system of prolate molecules (mesogens) sandwiched between flat 
walls, and exposed to a time-dependent external field. We find that a switch- 
able smectic domain forms for sufficiently high frequency. Although pressure 
and temperature are too low to induce an equilibrium smectic phase, the fluctu- 
ations of the pressure in the NEQ steady state match the pressure fluctuations 
characteristic of the equilibrium smectic phase. Furthermore, the wall-normal 
pressure fluctuations give rise to a tangential "fluctuation-vorticity" tensor that 
specifies the symmetry-breaking direction of the smectic layers. Our calculations 
demonstrate a novel method through which nanomaterials with a high degree 
of molecular order may be manufactured in principle. 

In nature, systems out of equilibrium are the rule and not the exception jg]. Yet, only 
recently, with the growing interest in nano- and mesoscopic phenomena, NEQ thermody- 
namics has attracted a remarkably growing interest jl, 9]. Some profound and pioneering 
results are already available 10Ml3|. but a coherent picture of NEQ physics is still lacking. 

A defining difference between the physics of NEQ and equilibrium systems is the presence 
of nonvanishing currents in the former, which are maintained by mechanical, thermal or 
chemical driving forces 2], [ljj]. In fact, many works 
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19] have addressed the statistical 



properties of these currents. Nonetheless, their physical role in a fluid is still not clear. 

Here, we perform molecular dynamics simulations of N = 4000 Gay-Berne-Kihara (GBK) 
molecules in presence of two atomically smooth flat walls. The GBK model has been success- 
fully used to reproduce a typical equilibrium liquid crystal (LC) phase diagram for prolate 
mesogens 20[. The isobaric-isothermal ensemble was used to avoid unphysical stresses on 
the system due to the combination of a cubic simulation box and the formation of anisotropic 
phases. 

The interaction of the molecules with the walls is modeled with a Lennard- Jones potential. 
Because of the presence of the walls the pressure in the system depends on the direction, 
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that is, it is a second-rank tensor, and not a scalar 

1 N 

V=-J2 [ mv * ® Vi + r 4 ® f, (1) 

i=l 

where Vj is the velocity of the z'th particle, m its mass, r, its position vector, fj the total force 
acting on it, V the total volume of the system, and the operator cg> represents the dyadic 
product. Further, because of the planar geometry of the system, only the diagonal Cartesian 
components V xx , V yy and V zz are nonzero, and the hydrostatic pressure is P = (V xx ) = (V yy ) 
(see Methods). 

In the case of anisotropic molecules, it is important to specify their preferential align- 
ment at the walls. A suitable quantity is the so-called "anchoring function" g(u) which 



discriminates energetically the orientation u of a LC molecule with resp ect to a surface 2 If . 



effectively defining a preferential direction, also called "easy axis" [22]. We consider a sys- 
tem confined by walls whose surface properties change periodically with time. Specifically, 
we give a temporal dependence to the anchoring function g(u, t) = A [u x V x + u y V y + u z V z ), 
where V x = V y = sin(a;t), V z = cos(ut), A is a constant, u is the angular frequency of the 
sinusoidal external field, and t is time. The effect of the external field is then to rotate the 
walls' easy axes with time. The two easy axes rotate in phase. Time- dependent, responsive 



surfaces are a growing field of research 23] . Clare et al. 24j have demonstrated that any 
specific anchoring of LCs can be selectively obtained by grafting semifluorinated organosi- 
lanes onto a surface. Further, the realization of time varying decorated surfaces has been 
reported where variable pH 25, [26], electrochemical 



properties 



or UV-light irradiation 



are used gents to effect the time dependence 

We present results for simulations in a range of temperature T = 4.0 — 6.0 at fixed 
P = 3.6 and investigate the behavior of the system as the period of the external field 
r = 2ir/u is varied. We use dimensionless units throughout this work (see Methods). We 
first analyze the dependence of the energy fluctuations of the fluid on r. Figure [1] shows 
that ((Auff) 2 ) grows linearly with r at large r (which is the usual behavior of the energy 
fluctuations for a sinusoidal field). At r 7 x 10 4 there is a crossover to a plateau. Because 
((Auff) 2 ) is proportional to the specific heat of the fluid, this crossover reflects a structural 
rearrangement of the system occurring at small periods (high frequency). 

Visual inspection of the molecular configurations reveals (Fig. [2^) that the molecules in 
the central portion of the fluid assemble in a well defined smectic domain (SD). This domain 
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FIG. 1: Crossover of the energy fluctuations upon formation of a NESS. Dependence of 
the fluid-fluid energy fluctuations on the period of the external field r. The dashed lines represent 
linear fits for large and small values of r. 

develops after a relatively small number of cycles of the external field and persists as long as 
the field is switched on; the SD disappears quickly once the external field is switched off. The 
fluid becomes heterogeneous. The molecules directly in contact with the walls rotate as the 
anchoring associated with the field changes. Between the contact layer at the walls and the 
SD there is a very turbulent layer, which does not show any spatial nor orientational order. 
Close observation of the molecular dynamics shows that individual molecules constantly 
leave this turbulent layer to join the SD or vice versa; however, the SD remains a stable 
feature. This is an instance of a nonequilibrium steady state (NESS). Only the molecules 
in the contact layer (i.e, the layer closest to the wall) are subject to direct interaction with 
the wall because of the short-range nature of the fluid-substrate potential U f a (see Fig. |2b 
and Methods). Therefore, the formation of SD must not be confused with a confinement 
effect, because the SD forms in the region where Uf s ~ 0. Hence, we believe the formation 
of the SD to be a general consequence of a time-dependent external field, irrespective of the 
precise realization of this field, so that our results are relevant to a broad class of physical 
situations. 

Now, a natural question to ask is: Why should a SD form? In equilibrium, the lowest P 
for which a smectic phase forms at T = 4.0 is more than twice the value of P investigated in 
this work. Thus, the equilibrium phase transition is too far removed to play any role here. 
It is also important to note that even the local value of the pressure in the region where 




FIG. 2: Molecular configuration of the simulated fluid and comparison with the exten- 
sion of the fluid-substrate interaction, a, Snapshot showing the lateral view of the system in 
the xz-plane at P = 3.6, T = 4.0 and r = 3 x 10 4 . The molecules closest to the walls (top and 
bottom layers) rotate with the external field. A smectic domain is clearly visible in the center of the 
system, b, Plot of the fluid-substrate interaction Uf s (see Methods) showing that it is effectively 
different from zero only in a region less than 5 molecular diameter in size. 

the SD forms is too low to explain a smectic state. The inset in Fig. [3] shows that the local 
pressure integrated over the SD volume, V zz , is too low compared with the same quantity 
but calculated for an equilibrium smectic state at the same P and T. 

We then consider the temporal fluctuations of the pressure ([AV zz (z)] 2 ) = ([P zz (z, t)] 2 ) — 
([Pzz( z , t)}) 2 . Figure [3] shows the dependence of ([AV zz (z)} 2 ) on the position along the z- 
axis. The confining walls are located at z/L z = ±1/2, where the pressure fluctuations are 
very large due to the molecular rotation. As we move towards the center of the system, 
([AV zz (z)} 2 ) decreases rapidly until it reaches an almost constant value. It is the main 
observation of this study that when a SD forms the NEQ pressure fluctuations match the 
value of the equilibrium pressure fluctuations in a smectic phase 

([AV zz (z)] 2 ) NEQ = ([AV zz (z)] 2 ) EQ . (2) 

Also, the region of the plateau of ([AV zz (z)} 2 ) coincides with the location of the SD. As 
t increases the SD shrinks and becomes less coherent; this correlates very well with the 
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FIG. 3: Normal pressure fluctuations. Dependence of the pressure tensor fluctuations on 
position along the z-axis for NEQ simulations at P = 3.6, T = 4.0 for r = 3 x 10 (0)> t = 7 x 10 4 
(□), and r = 9 x 10 4 (O). The black line shows the value of the pressure tensor fluctuations for 
an equilibrium simulation in the smectic phase. Inset shows the local pressure integrated over 

the SD volume V zz ( ). The NEQ value oscillates with time following the external field, and, 

interestingly, its average value is close to the value characteristic of an equilibrium isotropic phase 

( ), which is the thermodynamic equilibrium state at this P and T. The equilibrium value 

characteristic of the smectic phase is much larger (-•-•). From this we conclude that the local 
pressure is not large enough to drive the formation of a SD. 

behavior of ([AV zz (z)} 2 ) in Fig. [3] for r ^ 6 x 10 4 . This value deviates increasingly from 
([AV zz (z)] 2 ) EQ . 

The instantaneous value of V zz (z,t) may be treated as a stochastic variable resulting 
from the chaotic molecular motion and the oscillatory behavior at the walls. Therefore, 
to rationalize the coincidence between the SD formation and equation fl2]) we turn to a 
statistical description of the pressure profile in terms of a Fokker-Planck equation for the 
probability I%t) = (5(V zz (z,t) - p)) 

= -^[C( P )Uip,t)} + ±*L[D(p)ILM] (3) 

where 5(x) is the Dirac 5-function, (■ • •) represents the average over the molecular noise, 
Civ) = (Ap)/At = (AV zz )/At, D{z) = ({Ap] 2 )/At = ([AV zz ] 2 )/At in the limit At 



29(. Now, it is readily seen that C(p) vanishes because the field configuration is symmetric 



and therefore the transition probability is symmetric in the increment AV ZZ . Hence, the 




FIG. 4: Tangential pressure fluctuations. Dependence of the pressure tensor fluctuations 
([AVxxiz)] 2 } (black) and ([AT'yy(z)] 2 ) (red) on position along the z-axis for NEQ simulations at 
P = 3.6, T = 4.0 and r = 3 x 10 for a SD configuration with the layer-normal parallel to the 
x-axis (a), and at 45° with the x-axis (b). Insets show the top view of the cross sections of the 
two systems with the SD corresponding to the tangential pressure fluctuations calculated here. 

probability II (p, t) is governed not by the value of V zz but rather by the fluctuations similar 
to ordinary Brownian motion. 

Because the SD is a NESS we assume that mechanical stability is locally valid in the 
central part of the fluid, sufficiently removed from the walls. Locally then, 

V ■ V = 0. (4) 

In equilibrium, from equation fl3j) follows that V zz (z) = const, in the entire fluid. From the 
fact that ([AV zz (z)] 2 ) does not depend on z in the central portion of the fluid (where the 
SD forms, see Fig. [3]) we are led to assume that a similar relation to equation flj]) is valid 
for the pressure fluctuations 

V • (AV) = 0. (5) 

Similar to standard hydrodynamics, we can then define a "fluctuation-vorticity" tensor 
associated to the pressure fluctuations, u) = V x A7 3 . Because of the planar geometry 
of our system the pressure tensor components are only functions of z, such that u xx = 
— d(AVyy) / dz and u yy = d(AV xx )/dz. A larger slope of AV XX or AV yy then implies a larger 
ujyy or u xx , respectively. To test whether uj has physical significance we consider two systems 
with the same normal pressure fluctuations, i.e. the same r, but with different smectic-layer 
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normals. Figure shows the pressure-fluctuation profile {[AP aa (z)] 2 ), a = x,y, for a 
system exhibiting a SD with a smectic-layer normal parallel to the x-axis. The SD extends 
in the region \z\/L z < 0.2. In the same region AV XX has zero slope, while AP 9!/ exhibits a 
large slope. This, in turn, implies a vanishing u yy and a large u xx . The relative magnitude 
of u xx and u yy correlates with the orientation of the smectic-layer normal. Further, in Fig. 
[4b we show the pressure fluctuation profile for the second SD whose smectic-layer normal 
is at an angle of 45° with the x-axis. The two curves coincide (within numerical accuracy) 
indicating equal tangential components of uj. Therefore, from Fig. [4] we conclude that the 
fluctuation- vorticity uj determines the symmetry breaking direction of alignment of the SD. 

To conclude, we find evidence from NEQ computer simulations that pressure fluctuations 
can be easily tuned to drive a fluid system to a far-from-equilibrium state. The role of current 



fluctuations has been recognized 



17| as a stochastic variable characterizing NESS. Here, 



the physical picture emerging is that fluctuations in the momentum current (i) determine 
the NESS, and (ii) give rise to a secondary field that breaks the rotational symmetry in the 
x|/-plane. 

Self-assembly of molecules or supramolecular particles into layers, membranes, and vesi- 
cles is revolutionizing our control of matter across multiple length scales with far-reaching 



tuned to obtain the desired features 



applications in nanofluidic devices 22 . l30l . l3l| . Chemico-physical properties are carefully 



321 ] . However, in most cases they do not have any 



temporal dependence. The richness of NEQ phenomena in simple systems may suggest that 
combining the powerful new techniques of nanoscopic control with the application of time 
dependent external fields (temperature, electric or magnetic field, pressure and pH) may 
usher new ways to induce molecular self-assembly and even to simplify known tasks. In 
particular, the vorticity field may be used in the future to control the orientation of the 
ordered smectic domains which could be useful to manufacture new nanoscopic materials 
with distinct materials properties. 

Financial support from the International Graduate Research Training Group 1524 is 
gratefully acknowledged. 
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METHODS 



The fluid-substrate interaction is modeled with an "integrated" Lennard- Jones potential 
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(6) 



where ej s = 1 and p s a 2 = 2 2 / 3 tt is the areal density of a single layer of atoms arranged 
according to the (100) plane of a face-centered cubic lattice. The diameter a of the substrate 
atoms is equal to the LC molecular diameter. The quantity is the minimum distance 



331 ] between a LC molecule and the substrate located at z = —L z /2 (k = 1) and z = +L z /2 



(k = 2). The time-dependent anchoring g(u, t) is included in the attractive part of the fluid- 
substrate interaction. Dimensionless units are used throughout, that is, length is expressed 
in units of a, temperature in units of e^/kB, time in units of (cr 2 m/eg) 1 / 2 using m — 1, and 
pressure P in units of cr 3 /eg, where eg is the fluid-fluid interaction energy scale of the GBK 
model Q. 



We use a velocity- Verlet algorithm for linear molecules 



34j . and the simulations are car- 



ried out in the N 
Hoover barostat 



'T ensemble using a Nose-Hoover thermostat 35|, [36| and an anisotropic 



371 ]. whereby L z is kept fixed, while L x and L y are allowed to vary inde- 



pendently from each other, resulting in equal lateral average values of the pressure tensor 

(Pxx) — \Pyy)- — 

We use the "method of planes" 38] to compute the component V zz of the pressure tensor. 
Unfortunate this method cannot provide V„. and V„ by construction, which .ustead we 
compute following Harasima's method |39| . 
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